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Abstract 

Variations of frequency separations of low-degree p-modes are studied over the 
solar activity cycle. The separations studied are obtained from the frequencies of low- 
degree p-modes of the Global Oscillation Network Group (GONG). 10.7 cm radio flux 
is used as an index of solar activity. Small separations of the p-mode frequencies are 
considered to be mainly dependent on the conditions in stellar interiors. Thus they 
could be applied to diagnose the changes in the stellar interior. Our calculation results 
show that the magnitudes of variations of the mean large separations are less than 
1 a over the solar activity cycle. Small separations show different behaviors in the 
ascending and descending phases of activity. In the ascending phase, variations of the 
small separations are less than 1 a . However, the small separations have systematic 
shifts during 2004 - 2007. The shifts are roughly 1 o or more. The variations of the 
ratios of the small to large separations with time are similar to the changes of the 
small separations. The effects of the changes in the large separations on the ratios are 
negligible. The variations of the separations may be a consequence of the influence 
from the surface activity or systematic errors in measurements or some processes 
taking place in the solar interior. 
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1. Introduction 

It is well-known that the solar p-mode frequencies vary with the solar cycle (Woodard 
& Noyes 1985; Elsworth et al. 1990; Libbrecht & Woodard 1990; Elsworth et al. 1994; Jimenez- 
Reyes et al. 1998; Howe et al. 1999; Chaplin et al. 2007). It has been shown that the variations of 
the frequencies are significantly related to magnetic activities near the solar surface (Woodard 
et al. 1991; Jimenez- Reyes et al. 2004; Toutain & Kosovichev 2005). However, a number of 
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attempts have been made to detect the relation between the frequency shifts and the variation 
of the solar interior. It has been proposed that the frequency shifts may relate to changes 
taking place in the base of the convection zone (Jimenez- Reyes et al. 1998; Serebryanskiy & 
Chou 2005), or in the deep core (Jimenez- Reyes et al. 1998). 

Low-degree p-modes can be used to study the solar core because they can penetrate the 
deep interior (Dziembowski & Goode 1997). However, even these modes are primarily sensitive 
to the solar envelope structure. Thus frequency separations are proposed for a diagnostic 
tool. The frequency separations of low-degree p-modes in solar-like stars often proposed for 
diagnostic purposes are large separations A;(n) and small separations du + 2{n). The large 
separations approximate to the characteristic frequency z/o of a star, defined by 

Ai(n) = v ni i - v n _ u ~ v , (1) 



where 



u = (2 i R dr/c)-\ (2) 



o 



in which c is sound speed at radius r and R is some fiducial radius of the star. The small 
separations are sensitive to the structure of stellar core (Gough & Novotny 1990), defined by 

du+2{ n ) = v n,l — z/ n-l,/+2 /„\ 

- — ~ — VO, 

where ifi is the internal phase shift. For the low-degree p-modes, the <pi is mainly dependent on 
conditions in the stellar core and is almost independent of those outside the core (Roxburgh & 
Vorontsov 2003). These frequency separations have been extensively investigated by many au- 
thors (Christensen-Dalsgaard 1984; Ulrich 1986; Ulrich 1988; Gough 1987; Gough 1990; Gough 
2003; Gough & Novotny 1990; Roxburgh & Vorontsov 2003; Audard & Provost 1994; Roxburgh 
2005; Oti Floranes et al. 2005). Additionally Roxburgh (1993) defines differences doi(n) as 

(4) 



doi(n) = 



and Yang & Bi (2007) define differences cj-u+i (n) of low-degree p-modes as 

(Tl-ll+l(n) = -V n ,l-l + 2Vn,l-Vn,l+l /.n 
(p l+1 -\-<pi- 1 -2ipi \ ) 

7T u 

The differences 002 (w) are similar to the scaled small separations G^ +2 (n)/(2Z-|-3) in some cases 
and are mainly sensitive to conditions in stellar core. The difference (X02 averaged over n is, 
however, more sensitive to changes in the central hydrogen abundance of solar-like stars than 
the scaled small separations (Yang & Bi 2007). 

Moreover, the ratios dj/ +2 /Aj are considered to be essentially independent of the structure 
of the outer layers of a star and determined by interior structures (Ulrich 1986; Roxburgh & 
Vorontsov 2003; Roxburgh 2005; Oti Floranes et al. 2005). The ratios <i i/A and <7 2/A also 
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mainly rely on the interior structures (Roxburgh & Vorontsov 2003; Yang & Bi 2007). Thus 
changes in stellar interior could be indicated by the variations in these separation ratios. 

However frequency variations with solar activity are relative to the angular (I) degree, 
the azimuthal (m) degree, and the order n. The m-dependence is attributed to the solar 
near-surface activity (Jimenez- Reyes et al. 2004; Toutain & Kosovichev 2005) and could be 
removed by using the frequency centroids which incorporate all the m components of a mode 
(Chaplin et al. 2005). For the low-degree p-modes, centroid shifts can be regarded as being 
/-independent (Chaplin et al. 2005). Thus separations of frequency centroids of low-degree p- 
modes should be insensitive to effects of surface activity. Using the p-modes of the Birmingham 
Solar-Oscillations Network (BiSON) collecting data using a Sun-as-a-star technique, Chaplin 
et al. (2005) studied the impact of the solar activity cycle on the frequency separation ratios 
and found that the ratios change with the shifting level of global solar activity. Chaplin et al. 
(2005) pointed out that some m components are suppressed in the Sun-as-a-star observations 
and hence the ratios are expected to be affected by the surface activity. In this paper, we 
investigate the variations of frequency separations of low-degree p-modes as solar activity cycle 
proceeds using data collected by the resolved-Sun observations. In Section 2 we present our 
data and results. Then discussion and conclusion are represented in Section 3. 

2. Data and results 

In this work, we used the solar p-mode frequencies of the GONG (Harvey et al. 1996) 
to study the separations. The GONG was planned for measurements of medium and high 
spherical degrees, but low-degree p-modes can be measured too. Gavryusev & Gavryuseva 
(1999) showed that the low-1 p-modes of the GONG are in good agreement with that of the 
Michelson Doppler Imager (MDI) and the Global Oscillations at Low Frequency (GOLF). The 
raw data of GONG are measured in a time series of 36 days labeled a GONG month. A 108-day 
time series is constructed by concatenating 3 consecutive GONG-month time series. GONG 
Month Mode Frequencies (GMFs) were then estimated from the 108-day spectra (Anderson et 
al. 1990; Hill et al. 1998; Komm et al. 1999; Howe et al. 1999). In the observation of GONG, 
the solar disk is imaged onto many pixels. The obtained images can then be decomposed into 
their constituent spherical harmonics. For each degree / and overtone n, this imaging strategy 
provides access to all 21 + 1 components with each being tagged by an azimuthal degree m 
(Chaplin et al. 2004). 

The frequencies, ^ n ,«,m, extracted from resolved-Sun observations are usually represented 

by 

21 

v n ,i,m = < jZ + ^a i (n,/)/P/(m), (6) 
where V\[m) are polynomials related to Clebsch-Gordan coefficients (Ritzwoller & Lavely 1991; 
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Fig. 1. Differences between the averages obtained using equation (7) and the centroids estimated by the 
GONG team. Circles represent the differences of frequencies of the GONG months 8, 9, and 10 (low 
activity), while crosses correspond to the differences of frequencies of the GONG months 58, 59, and 60 
(high activity). Triangles show the errors of GONG centroids of month 8. 

Chaplin et al. 2004) and v° l is the so-called central frequency of the multiplet (i.e. frequency 
centroid). Using the polynomials V\{m) given by Chaplin et al. (2004) and equation (6), one 
can get the central frequency of low-degree p-modes of the resolved-Sun observations, 

m=—l 

The central frequencies are sensitive to the spherically symmetric component of the internal 
structure of a star and are usually used as inputs to inversions of the internal structure (Chaplin 
et al. 2004; Chaplin et al. 2005). 

In order to study the variations of frequency separations with the solar activity cycle, 
we used the GMFs observed between 1996 and 2007. Centroids v n ,i and individual frequencies 
v n ,Lm of the GMFs were published by the GONG team 1 . Some centroids u n> i published are 
inconsecutive in the order n, especially for the modes with 1 = 1, but the frequencies v n ,i,m are 
consecutive. In order to obtain separations, we used the averages v° z of all 2/ + 1 components 
of each mode instead of the centroids z/ n ,i- ^ n fig ure 1? we compared the averages obtained 
using equation (7) with the available centroids estimated by the GONG team. The differences 
between the frequencies obtained using equation (7) and the available GONG centroids are 
almost much less than the errors of the GONG centroids in the range of 2400 /xHz < v < 3500 
//Hz except a few modes, whose differences are close to the errors of the GONG 
centroids. In this range, the averages obtained using equation (7) are good consistent with the 
centroids estimated by the GONG team. All the following frequency separations are obtained 
from the averages. 

Figure 2 shows fractional changes of the separation ratios in two sets of GMFs, one at 

1 ftp://gong2.nso.edu/MFS/ 
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Fig. 2. Fractional changes between the separations ratios in two GMFs, one at a high level (2001) and 
one at a low level of activity (1996) (the differences are in the sense "high data - low data"). 

Table 1. Results of linear fits between separations and time. Left panel for the ascending phase of so- 



activity, while 


right panel for the 


during < 




2004 to 


2007 except for the 


large separations. 


Ascending phase 


From - 2004 to 2007 


Separation 


Slope (nHz/year) a 


x 2 




R 


Slope (nHz/year) a 


x 2 


R 


< d Q2 > 


- 6.9 ± 5.1 


0.228 




0.177 


- 35.4 ± 7.6 


0.019 


0.697 


< d 01 > 


1.1 ± 5.3 


0.244 




0.027 


- 76.9 ± 13.4 


0.058 


0.768 


< C 2 > 


- 9.9 ± 7.1 


0.446 




0.184 


108 ± 23 


0.179 


0.692 


< di3 > 


- 8.6 ± 3.3 


0.092 




0.327 


43.0 ± 9.0 


0.029 


0.691 


< A > 


11.8 ± 2.8 


0.072 




0.485 


- 14.7 ± 2.3 b 


0.041 


0.654 


< Ai > 


5.2 ± 2.2 


0.041 




0.303 


- 11.4 ± 1.3 b 


0.012 


0.771 


< d 02 /3A > 


- (1.9 ± 1.2)xl0" 5 


1.3x10" 


-6 


0.205 


- (8.8 ± 1.9)xl0~ 5 


1.2xl0~ 7 


0.698 


< doi/Ao > 


(5.6 ± 3.9)xl0~ 5 


1.3x10" 


-5 


0.019 


- (5.7 ± 1.0)xl0- 4 


3.2xl0~ 6 


0.766 


< (T02/A0 > 


- (7.5 ± 5.3)xl0~ 6 


2.5x10" 


-5 


0.188 


(8.0 ± 1.7)xl0- 4 


9.8 xl0~ 6 


0.694 


<di 3 /5Ai > 


- (1.4 ± 0.5)xl0~ 5 


2.0x10" 


-7 


0.348 


(6.5 ± 1.4) xlO -5 


6.2 xl0~ 8 


0.699 



a. Not for ratios. 



b. The results were obtained from the data between ~ 2001 and 2007. 

a low level and one at a high level of activity. The error bars indicate 1 a errors, which were 
propagated from uncertainties on the individual frequencies of the GMFs. The departure of 
the separation ratios from the null level is within 1 a but shows a systematic shift, which is 
consistent with the prediction of the 108 day spectra of Chaplin et al. 2005 (in their figure 1). 
However, just as the observed shifts of ratios of the BiSON frequencies (Chaplin et al. 2005), 
'the shifts are clearly not significant for individual estimates of the ratios'. In the interest 
of temporal variations of separations, for each set of separations with the same degree /, we 
averaged the separations over n in the range of 2400 /iHz < u nj i < 3500 /iHz, which was implied 
by symbols <>. The frequency range is often used to obtain a mean frequency shift (Elsworth 
et al. 1994; Jimenez- Reyes et al. 2004). 

Figure 3 shows the temporal variation of mean separations < du + 2 >, < 002 >, < c?oi >j 
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Fig. 3. Temporal variation of separations. The error bars represent 1 a errors. The thin lines show the 
linear fits between the plotted quantities in the ascending phase of the solar activity, while the thick lines 
indicate the linear fits in the descending phase of the solar activity. Sunspot number and 10.7 cm radio 
flux are plotted too. 1 SFU = 1CT 22 W • m 2 • Hz" 1 . 
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Table 2. Results of linear fitting between separations and 10.7 cm radio flux. 



Ascending phase 


From - 2004 to 2007 


Separation 


Slope (nHz/SFU) 


x 2 


R 


Slope (nHz/SFU) 


x 2 


R 


< d 02 > 


- 0.34 ± 0.19 


0.223 


0.235 


2.13 ± 0.69 


0.026 


0.542 


<d i > 


0.024 ± 0.197 


0.244 


0.016 


5.6 ± 1.1 


0.068 


0.722 


< C02 > 


- 0.33 ± 0.27 


0.449 


0.162 


- 8.4 ± 1.8 


0.176 


0.699 


<d 13 > 


- 0.25 ± 0.13 


0.096 


0.254 


- 3.2 ± 0.8 


0.032 


0.657 


< A > 


0.33 ± 0.11 


0.081 


0.361 


0.57 ± 0.09 a 


0.041 


0.648 


< Ai > 


0.22 ± 0.08 


0.040 


0.329 


0.45 ± 0.05° 


0.012 


0.791 



a. The results were obtained from the data between ~ 2001 and 2007. 

< Ao > and < Ai >. Sunspot numbers and 10.7 cm radio flux 2 (F10), which are always used 
as an index of the solar activity, are also plotted as a function of time. The activity indices 
have been averaged over the same time interval covered by the GONG date. The error bars 
show 1 a errors, which come from the formal errors on the frequencies. When 
averaging these separations from each spectrum, the correlations of neighbouring 
separations were considered. The lines in this figure represent the results of linear least- 
squares fits. The fits in the descending phase of solar activity were divided into two periods 
in accordance with the variations of separations: one from ~ 2001 to 2004 and one from ~ 
2004 to 2007. Variations in the latter period are more obvious and systematic. The mean 
large separations < Ao > and < Ai > increase/decrease in the ascending/descending phase of 
solar activity and reach a maximum in 2001. However, the magnitude of the shifts of the large 
separations is about 0.07 /iHz, which is less than 1 a. Variations of small separations with 
time are more complicated than those of the large separations. In the ascending phase of solar 
activity, the separations < du +2 >, < cr 2 > and < rf i > almost have no systematic shifts; the 
magnitudes of the shifts of these separations, obtained from their regression equations, are also 
less than 1 a. In the descending phase of solar activity, variations of these separations seem to 
be disorderly. Between 2001 and 2002, the separations < 002 > and < d\% > increase suddenly. 
They decrease firstly and then increase obviously between ~ 2001 and 2007. However, the 
magnitudes of their shifts are roughly 1 d or more. The separation < d m > has an obvious and 
systematic decrease between ~ 2004 and 2007. The magnitude of the variation is about 1.5 a. 
The separation < d 02 > decreases slightly between ~ 2001 and 2007, but the magnitude is less 
than 1 a. The results of linear fits: slope, residual sum of squares x 2 ; and correlation coefficient 
R of regression equations, are summarized in Table 1. 

In order to further understand the variations of separations with the solar activity, 
Figure 4 shows separations plotted against F10 in the ascending phase of solar activity. Figure 

2 www.ngdc.noaa.gov/stp/SOLAR/getdata.html 
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no (s.f.u) no (s.f.u) 

Fig. 4. Mean separations are plotted as a function of 10.7 cm radio flux in the ascending phase of solar 
activity. The error bars indicate 1 a errors. The lines represent the best linear fits between the plotted 
quantities. 

4 shows that the separations almost not vary with activity. The results of linear least-squares fits 
between the separations and F10 are summarized in the left panel of Table 2. The magnitudes 
of variations of these separations are only about 0.04 /iHz, which is less than 1 a. Moreover, 
Figure 5 shows that the changes of the large separations with the F10 in the descending phase 
of the solar activity. The shifts of the large separations are also less than 1 a. Figure 5 also 
shows the variations of small separations with the F10 during ~ 2004 - 2007. In this period of 
time, the changes of these separations with the F10 are obvious and systematic. However, the 
magnitudes of variations of these separations are around 1 a or more. The results of linear fits 
between the separations and F10 are shown in the right panel of Table 2. 

In Figure 6, we show the temporal variations of separation ratios < <io2/3Ao >, 
< 6^3/5 >, < doi/A > and < 002/ A >. The results of linear fits between the ratios and 
time are presented in Table 1. The changes of the ratios with time are similar to those of the 
corresponding small separations. In the ascending phase of solar activity, shifts of the ratios 
are less than 1 a. In the descending phase of activity, from ~ 2004 to 2007, the changes of the 
ratios with time are systematic. However, the magnitudes of the changes are still roughly 1 a 
or more. Since the fractional changes of the large separations are much less than those of small 
separations during the activity, the contribution to the changes in the ratios from the variations 
of the large separations is much less than that from the small separations. The effects of the 
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Fig. 5. Mean separations observed between ~ 2004 and 2007 are plotted as a function of 10.7 cm radio 
flux. But the large separations are observed between ~ 2001 and 2007. The error bars indicate 1 a errors. 
The lines represent the best linear fits between the plotted quantities. 




Fig. 6. Temporal variations of separation ratios. The error bars indicate 1 a errors. The lines show the 
linear fits between the plotted quantities. 
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changes in the large separations on the changes in the ratios are almost negligible. The changes 
in the ratios with time are mainly determined by the changes in the small separations. 

The magnitudes of temporal variations for some individual separations are very large 
but should not be significant. In the ascending phase of the solar activity, the magnitudes of 
variations of separations, obtained from regression equations, are about 0.04 //Hz, which is less 
than 1 a. During ~ 2004 - 2007, the changes in separations are systematic. However, the 
magnitudes of the changes are roughly 1 a or more. Using BiSON data, Chaplin et al. (2005) 
show that variations of the separation ratios between low and high activity is at a marginal 
level of significance. Moreover, using simple averages, from the gradient of mean frequency 
shift of the MDI given by Jimenez- Reyes et al. (2004) in their Table 1, one can get a shift 
gradient in the centroids would be 1.61 and 2.02 nHz/SFU for / = and 2 respectively. For a 
120 SFU increase in solar activity, the change in separation < d 2 > would then be about 0.05 
/iHz, which is consistent with our result. 

3. Discussion and conclusion 

The characteristic frequency u Q , i.e. the large separation, is an integral property of 
the entire star. Because the solar sound speed is a decreasing function of the solar radius, the 
contribution to v$ from the outer layers of the Sun is larger than that from the solar inner regions 
(Gough & Novotny 1990). Thus the large separation should be more sensitive to the changes 
in the outer layers than those in the inner regions and could be affected by the near surface 
activity. Moreover small separations are more sensitive to the conditions in stellar interior than 
those in outer layers. Thus the large separations should be more sensitive to the activity than 
the small separations. However, in fact, the magnitudes of variations of the large separations 
are about 0.07 /iHz during the activity, which is less than 1 a. The magnitude of variations of 

< do2 > is also less than 1 a over the activity cycle. However, other small separations show a 
systematic shift during ~ 2004 - 2007, and the magnitudes of their variations are roughly 1 a 
or more. 

Small separations < du + 2 >, < O02 > and < cfoi > show a sudden change between 2001 
and 2002. The fractional changes shown in Figure 2 of ratios in two GMFs (one in 1996 and 
one in 2001) are consistent with the prediction of Chaplin et al. (2005). However, if we use the 
data of 2002 instead of the data of 2001, the results would be different. The small separations 
show different behaviors in the ascending and descending phases of the cycle. In the ascending 
phase, variations of the small separations are less than 1 a. In the descending phase, the small 
separations show a systematic shift during 2004 - 2007. And the magnitudes of the changes in 

< di3 >, < O02 > and < doi > are roughly 1 a or more. It is not clear whether the variations 
are caused by the systematic errors in measurements. 

The large separations increase slightly with increasing of the solar activity. But the 
magnitudes of temporal variations of the mean large separations of the low-degree p-modes are 

10 



only about 0.07 /iHz over the activity cycle, which is less than 1 a. Small separations < du + 2 >, 
< (T02 > and < doi > show different behaviors in the ascending and descending phases of the 
cycle. In the ascending phase, variations of the small separations are less than 1 a. However, the 
small separations show a systematic shift during 2004 - 2007. The magnitudes of the changes 
in < doi >, < O02 > and < d%3 > in this period of time are about 1 a or more. The variations 
of the ratios of the small to large separations with time are similar to the changes of the small 
separations. Changes in the large separations have a negligible impact on the variations in 
the ratios. The variations of the separations may be a consequence of the influence from the 
surface activity or systematic errors in measurements or some processes taking place in the 
solar interior. For further work, the data over cycle 22/23 would be required. 
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